Characterizing and quantifying frustration in quantum many-body systems 
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We present a general scheme for the study of frustration in quantum systems. We introduce a 
universal measure of frustration for arbitrary quantum systems and we relate it to a class of entan- 
glement monotones via an exact inequality. If all the (pure) ground states of a given Hamiltonian 
saturate the inequality, then the system is said to be inequality saturating. We introduce sufficient 
conditions for a quantum spin system to be inequality saturating and confirm them with extensive 
numerical tests. These conditions provide a generalization to the quantum domain of the Toulouse 
criteria for classical frustration-free systems. The models satisfying these conditions can be reason- 
ably identified as geometrically unfrustrated and subject to frustration of purely quantum origin. 
Our results therefore establish a unified framework for studying the intertwining of geometric and 
quantum contributions to frustration. 

PACS numbers: 03.67.Mn, 03.65. Ud, 05.50.+q, 75.10.Jm 
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Frustration consists in the impossibility of determin- 
ing configurations that minimize simultaneously the en- 
ergy of competing interactions [H-Q. In recent years it 
has been realized that classically unfrustrated systems 
can have frustrated quantum counterparts 0d5|- Indeed, 
in the quantum case, additional sources of frustration 
may arise due to non-commutativity and entanglement 
[H, [§]. Hence, although the notion of frustration has of- 
ten been considered from the perspective of statistical 
physics, in the quantum domain interesting novel phe- 
nomena take place even in the presence of few entan- 
gled elements [B|, H, ■ It is however far from clear how 
to distinguish geometric from purely quantum frustra- 
tion and whether the distinction is free of ambiguities. 
To the best of our knowledge this issue has never been 
addressed. Moreover, despite frustration being a well- 
defined and intuitive concept, a measure for quantum 
systems which fully captures all its aspects is still lacking. 
Existing proposals for quantifying frustration in quantum 
systems focus on particular aspects of the phenomenon. 
Some proposals are based on classical equivalents of a 
given quantum system making it impossible to re- 
cover information about quantum correlations. Others 
apply only in the thermodynamic limit [ill [l2|, or on 
the competition between local and nonlocal interactions 
from a purely quantum perspective thus making it 
challenging cither to study finite-size systems or to un- 
derstand the transition to classically frustrated systems. 
The need for a systematic investigation of this issue at a 
foundational level is thus quite compelling. Even more 
so, as quantum simulators of classical as well as quantum 
frustrated magnetic systems, at least in the simplest Ising 
and J\-Ji cases, are being demonstrated with ultracold 



atoms in optical lattices [13!, LLJ] and trapped ions [15( , or 
proposed, e.g. with cold Coulomb crystals [3]. Further- 
more, such an investigation would be of great help and 
guidance for the design of entanglement-based numerical 



simulations of frustrated quantum spin models 17 1. 

The aim of this work is to introduce a universal mea- 
sure of frustration for quantum systems and define a uni- 
fied framework suitable to understand the intertwining 
of the geometric and quantum contributions to frustra- 
tion. To this end we focus on the microscopic properties 
of finite-size models from a purely quantum perspective 
without resorting to scmiclassical approximations. Af- 
ter introducing a universal measure of frustration, we 
prove that it is an upper bound to a class of entanglement 
monotones that in particular cases reduce to the bipar- 
tite geometric entanglement. We then establish sufficient 
conditions for a quantum spin system to saturate the 
bound and support them with extensive strong numeri- 
cal evidence. Finally, we show how these conditions es- 
sentially generalize to the quantum domain the Toulouse 
criterion for frustration-free classical systems. 

Consider a system in a pure state p = \^p){ip\, and let 
II be the projector onto a given subspace. Then 



/=l-tr[pll], 



(1) 



quantifies how much p fails to fully overlap with the sub- 
space selected by II. Let now p = \G)(G\ be the ground 
state (GS) of a many-body system H = J^s an d lis 
the projector onto the GS of the local interactions hs cor- 
responding to subsystem S. Then fs = 1 — trills] is a 
well-defined and unambiguous measure of the frustration 
of hs- On the other hand, denoting by R the rest of the 
system, consider the following entanglement monotone: 



E^(S\R) = l-J2^(Ps) 



(2) 



where Xf(ps) are the eigenvalues of ps = tr^p in decreas- 
ing order and d < dim[Hs], with Hs the Hilbert space of 
S. Notice that E^(S\R) = E^(R\S). The r.h.s. of Eq. 
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([2]) vanishes only on states with Schmidt rank smaller or 
equal than d [l8|; for d = 1 it reduces to the bipartite 
geometric entanglement, defined as the distance from the 
set of bi-separable states . Fix now d to be the de- 
generacy of the local interaction hs (i.e. II is rank-d). 
Then, by Cauchy interlacing theorem (See also Lemma 1 
in the supplementary material) it follows that 

fs > 4 d) . (3) 

where E^ = E^ d \S\R) is the distance from \G) to the 
closest state with Schmidt rank r < d [l8| . 

Despite its apparent simplicity, Eq. Q has the remark- 
able feature of directly relating frustration to entangle- 
ment. This quantitative relation holds for any pure state 
p and any interacting quantum system and hence is uni- 
versal. Actually, Eq. © holds as well for any mixed 
state, although in this case E s d ~* is no longer an entan- 
glement monotone. An immediate consequence of Eq. 
<j3j> is that the frustration-free condition, fs = OVS* is 
a bound on the maximum Schmidt rank of the global 
ground state \G). For interactions hs with nondegener- 
ate local GSs |Gs), d = lVS 1 , this implies the separabil- 
ity of the global GS in the tensor product of the local 
GSs: \G) = ® s \Gs)- That the absence of frustration 
should be related to some form of factorization of the GS 
had already been observed, at a semi-quantitative level, 
for models in transverse fields Q. On the other hand, 
saturation of the inequality, fs = Eg for some S, im- 
poses a block-diagonal form of the reduced state ps, with 
eigenvalues \\{p),i = 1, . . . , d corresponding to the block 
spanned by II (see Lemma 1 in the supplementary ma- 
terial). Hence ps cannot exhibit coherence between the 
lowest and excited energy levels of hs'- the largest contri- 
bution to ps must come from the local ground subspace. 

Summarizing the above discussion, we can state that a 
GS p = \G)(G\ of a many-body Hamiltonian of the form 
H = hs is a frustration-free (FF) state if and only if 
fs = VS 1 , and is an inequality saturating (INES) state 
if and only if fs = Eg VS. Clearly, a FF state is also an 
INES state. A Hamiltonian H is then a FF Hamiltonian 
if all its GSs are FF states, and is an INES Hamiltonian 
if all its GSs are INES states. In general, it is easy to 
show that states with at least one fs > have higher 
energy compared to the corresponding FF state. Hence, 
if a model is globally degenerate and one global GS is 
a FF state, then this is true for all other GSs. Unlike 
the FF property, the INES property is not universal, i.e. 
it does not necessarily apply to all the GSs of a system. 
This is due to the fact that unlike the FF property, the 
INES property does not specify the GS energy. Total 
frustration F can then be defined out of the set of local 
measures /s's by arithmetic or weighted averaging. 

We will now exploit Eq. ([3]) to generalize the classi- 
cal criteria for the absence or presence of frustration to 
the quantum domain and to understand the intertwin- 



ing of geometric and purely quantum contributions to 
frustration. According to Toulouse [Formulation 1] A 
classical Hamiltonian H is frustrated if and only if it 
is impossible to transform H into a fully ferromagnetic 
Hamiltonian by means of local spin inversions. This oc- 
curs only when a closed loop exists with an odd num- 
ber of antiferromagnetic interactions [H, [1, HJ. Indeed, 
the Toulouse criterion computes exactly the parity of 
antiferromagnetic bonds on a closed loop, according to: 
[Formulation 2] For a given Hamiltonian a loop is frus- 
trated if the quantity P = ( — l) Naf = —1, where N a f 
is the number of antiferromagnetic bonds. Next, con- 
sider the simplest quantum extension of classical models, 
the quantum Ising model H = J2z<j ^ij^i^j = 2). 
Clearly, the model is unfrustrated if the GS is of the form 

IV'> = ®l 1 (sfr>(a\n---) + m ■■■)), with 7i = 0,1. 

Indeed, in this case fij = E^ } = VS* = ij, and Eq. © 
is saturated by all GSs. On the contrary, if the model 
is frustrated there exists at least one GS which is not an 

(2) 

INES state. In fact, all separable GSs exhibit E-j = 
on all pairs and fij > on at least one pair. Therefore, in 
terms of Eq. ([3]), the immediate extension to the quantum 
Ising case of the Toulouse condition is [Formulation 3] : A 
quantum Ising Hamiltonian is frustrated if and only if it 
is not an INES model. Although the Ising model does not 
contain quantum features, such as non-commutativity of 
the local interactions, and therefore can only exhibit ge- 
ometric frustration, the possibility of restating Toulouse 
criterion in terms of / and an entanglement monotone 

(2) 

E s is remarkable as it provides the first bridge between 
the classical and the quantum domains. 

From the Ising example we learn that in the quantum 
domain the relevant information to detect frustration is 
the existence of GSs not being INES states, rather than 
that of GSs not being FF states. This is because for 
quantum systems frustration arises non only from topo- 
logical constraints (so-called geometric frustration) but 
as well from purely quantum ones. Consider the fol- 
lowing classically unfrustrated Heisenberg Hamiltonian 
for four spins on an open chain: H = S, • Sj+i, 
where S, = (Sf,Sf,Sf), d = 1 and the global GS 
of H is non degenerate. The quantum version of this 
model is frustrated, as in the GS the total measure 
of frustration F = N^ 1 J2ij fij = g(3 — Vo), where 
N = 3 is the number of bonds. Remarkably, the frus- 
trated GS is an INES state, hence the frustrated model 
is still an INES model, in contrast to what occurs in 
the Ising case. Let us now add to H some geometric 
frustration i.e. H' = H + J^. Sj • Si+2- Now one has 
fij = 2/5 > E\f = l/3V(ij), and hence the frustrated 
GS is not an INES state. Indeed, the Ising model can 
only exhibit geometric frustration, even in its quantum 
version, as the local GSs are always unentangled. On 
the contrary, the GS of each Heisenberg local pair term 
is a maximally entangled Bell state. Accordingly, the 
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model can either be an INES model or not, depending 
on the topology of the system. Actually, even if geomet- 
ric frustration is absent, yet monogamy of entanglement 
[iol ] prevents the global GS from minimizing all the lo- 
cal terms, and the model is both frustrated and INES. 
This strongly suggests that failure to saturate Eq. is 
a signature of the presence of geometric frustration. 

From the above discussion it follows that it would be 
highly desirable to identify conditions detecting a priori 
the nature of the frustration in a given quantum system. 
Considering for the moment being only models with non- 
degenerate global GS, we approach the problem by ob- 
serving that the classical Toulouse condition contains two 
main ingredients: a ferromagnetic model, which serves 
as the prototype of FF models, and a gauge group under 
which the FF property is preserved. We thus look for a 
prototype INES model and a gauge group under which 
the INES property is preserved. To fix the stage, consider 
the general XYZ (Heisenberg) exchange Hamiltonian 



H-- 



Ejfl qfi qfi 



(4) 



x, y, z, with coupling vectors Jjj = 



where \i 

{Jfj, Jfj, Jfj) on arbitrary graph geometries. 
[Prototype model] : A non-degenerate quantum Hamil- 
tonian as in Eq. Q will be called prototype if la) there 
exists at least one local GS common to all local pair in- 
teractions hij's, and lb) each local coupling vector «7y 
has non-negative components (i.e. it is a full ferromag- 
net: > Vi, j, fx and there exists a common two-body 
state vector \<j>) which is ground state of all hij). 
[Conjecture 1]: All prototype models are INES. 

Note the remarkable parallelism between condition la) 
and the FF condition, which essentially requires the exis- 
tence of at least one global GS common to all hij's. From 
our conjecture, it follows that if d — 1 then all h^s must 
admit the same Bell state as GS. Since the GS of a given 
hij with positive Jij is determined by the lowest value 
component of the vector, it follows that it must be along 
the same axis for all Jij in the prototype model. Any 
model obtained from the prototype by means of tensor- 
ing local unitary operations is clearly still an INES model 
with the same set of and E^^s. In fact, one has 
[Conjecture 2]: Any model H' = H Tk which can be ob- 
tained from a prototype model H by partial transposition 
on any set of sites K is still an INES model. 



Applying partial transposition [2l| to a prototype 



Hamiltonian changes and E^; however, saturation 
of Eq. §5§ is preserved. Although an analytical proof 
of these conjectures appears challenging, in the supple- 
mentary material we provide compelling numerical ev- 
idence of their validity. This is a very significant re- 
sult, as it implies that not only local rotations but also 
partial transpositions do preserve the INES property. 
It is not a priori obvious that this would be the case, 



since, as already noted, both sides of inequality Q arc 
typically changed by a partial transposition. Indeed, 
partial transposition is intimately related to parity, as 
Pk(H) = Sj[H Tk Sj{, where Pk is the parity transforma- 
tion on site k [P(S) = — S). Thus, according to Conjec- 
ture 2, the local gauge group is G = SU{2) ® Z 2 . An 
element of SU(2) acting on S^ is represented by a trans- 
formation R £ SO(3), whereas a parity transformation is 
simply — 1. Hence, spin operators transform according to 
the 0(3) representation of G, and a local gauge transfor- 
mation g = ^ i gi (with gi € G) maps two-body intcrac- 
tions E^ffi into E,J^ T ^rW Given 
a general XYZ Hamiltonian H = — J2ij^2/j, Jtj^i $j > 
we derive necessary and sufficient conditions for H to 
be equivalent to a prototype model under the action of 
some g S G® N . Consider two sites a and b and let 
p(a — >• b) — {(a, ii), (£2, £3), ■ ■ ■ , (£fc, £>)} be any path from 
a to 6, where all pairs (£, j) € p interact. Define the sign 
of the path p as 



n(p(a -^b)) = diag(s 



(5) 



(ij')6P 



where s^- = sign( J^) and we define sign(O) = 1, with the 
product taken over all adjacent pairs ij belonging to the 
path p(a —> b). Note that ir(p(a — > b)) — ir(p(b —> a)). 
We state the following 

[Theorem] Necessary and sufficient conditions for a 
Heisenberg Hamiltonian H to be mapped into a prototype 
model by a local gauge transformation are 

2a) All coupling vectors Jij, have the smallest absolute 
value component along the same axis; 

2b) For any pair of spins a and b, ir(p(a — > b)) is inde- 
pendent of the path p from a to b. 

As shown in the supplementary material, condition 2b) 
guarantees that the system can be brought to a fully 
ferromagnetic system. When this condition is met, the 
dependence of it on the path can be dropped. Then one 
can consider ir(a — > b) as a "conservative field" which 
dictates the local transformation that has to be applied 
in b so that the sign of any path from a to b is positive. 
If this holds for all b, one can turn all the couplings to 
positive by means of local transformations. On the other 
hand, condition 2a) is related to the existence of the same 
Bell state as the common ground state of all the two-body 
interactions. A remarkable simplification occurs when 
the system is translationally invariant. In that case: 
[Theorem] If a model satisfies ir(£) = +1 , for every 
loop £ in the elementary cell and all coupling vectors have 
the smallest component along the same axis, then it is an 
INES model (Proof in the supplementary material) . This 
Theorem encodes in compact form three Toulouse's crite- 
ria, one for each spatial direction. This strongly suggests 
that the class of INES models defined by conditions 2a- 
b ) can be identified as that of geometric frustration- free 
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spin-1/2 quantum models. Our analysis reveals that the 
quantum nature of the model affects the very notion of 
geometric frustration. Condition la) or equivalcntly 2a) 
is a consequence of the existence of three inequivalent fer- 
romagnetic states, i.e. the triplet states. This constraint, 
however, is only relevant in systems with inhomogeneous 
couplings. Otherwise, conditions 2a-b) simply reduce to 
a generalized form of Toulouse's criterion. Being the class 
of geometric frustration-free quantum systems more re- 
strictive than its classical analogue, one may expect that 
further investigation reveals generic properties which so 
far failed to be properly generalized. 

On the other hand, the class of geometric unfrus- 
trated quantum systems is strictly larger than that of 
frustration-free systems. A deeper investigation of ge- 
ometric FF quantum models might thus unveil several 
relevant applications. For example, being the largest 
eigenvalues of pij identified by ELy may provide guidelines 
for designing optimal entanglement renormalization algo- 



rithms and other tensor network ansatz 22j. Moreover, 
the measure of frustration that we have introduced and 
its relation to geometric entanglement is not restricted 
to local two-body Hamiltonians or spin-1/2 systems and 
might be exploited to gain understanding of geometric 
frustration in arbitrary quantum many-body systems. 

According to our results, the INES class of models is 
larger than that of geometric FF, thus suggesting that 
the INES property has deeper implications than merely 
detecting the presence or absence of geometrical frus- 
tration. For example, the fact that INES geometrically 
frustrated models may behave differently from non INES 
geometrically frustrated ones hints at a significant role 
of the INES property as a diagnostic tool for quantum 
phase transitions in complex and frustrated models. In- 
deed, consider a frustrated system such as the elemen- 
tary cell of the pyrochlore lattice H = J(Si • S2 + S3 • 
S 4 ) + (Si + S 2 ) ■ (S 3 + S 4 ) H. For J > 1 the sys- 
tem is INES since it is in a dimer phase and all dimers 
have / = E^ 1 ' = 0, whereas all other bonds are in a 
maximally mixed state, hence / = E^ = 3/4. Ex- 
act diagonalization shows that the transition from the 
dimer (J > 1) to the plaquette phase (0 < J < 1) cor- 
responds to the transition from INES to non- INES, with 



fi 



34 



1 > E 



E. 



(i) 
34 



2/3. This suggests 



a correspondence between the different quantum phases 
of a frustrated model and its INES or non-INES charac- 
ter. Indeed, the relation between geometric frustration 
and exotic matter phases has been previously pointed out 
One may also ask how the presence/absence of ge- 
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ometric frustration and the INES/non-INES nature of a 
system gets reflected in computational [24| , information- 
theoretic 13, 25 1 or thermodynamic terms [13]. Neces- 
sary steps toward further investigation of all these open 
fundamental questions will require achieving a rigorous 
mathematical control and understanding of our quan- 
tum Toulouse conditions as well as the generalization of 



our approach to globally degenerate systems. Although 
some instances of degeneracy can be easily accommo- 
dated within our scheme (e.g., odd number of spins), 
some others (e.g. thermodynamic degeneracy) are still 
elusive, but in very simple cases such as the XY model. 

The essential and intriguing role played by partial 
transposition cannot go unnoticed. The fact that the pro- 
totype models preserve their INES character under par- 
tial transposition is to be expected if INES has anything 
to do with the presence or absence of geometric frustra- 
tion. Nevertheless, the preservation of the INES property 
is far from trivial because, unlike with the case of local 
unitary transformations, the GSs of a Hamiltonian be- 
fore and after partial transposition are not immediately 
related. Indeed, we expect that the rigorous proof of our 
conjectures will shed further light on the role played by 
partial transposition and its relation to geometric frus- 
tration. From a directly physical point of view, as already 
noticed.our quantitative analysis confirms previous ev- 
idence p that some form of factorization of the global 
GS is a necessary ingredient in the characterization of 
FF systems. It is tempting to speculate that investigat- 
ing the intertwining between frustration and factorizabil- 
ity of higher order (fc-scparability, including dimcrization 
and trimcrization) could lead to a unified framework for 
the understanding of the relations between frustration, 
the role of hierarchical geometric entanglement [III] in 
collective quantum phenomena, and the characterization 
of entanglement and fc-separability by local unitaries [2(| . 
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Supplementary Material 

Lemma 1: Let II be a d- dimensional projector, II = 
t — II its complement and p a density operator. Then, 
saturation of the inequality tr[ j oII] < X^=i ^i(p) implies 
p = diag(pi,p 2 , ■ • -Pd) © IlpU. where pi = \\{p) are the 
eigenvalues of p in decreasing order. 

Proof: Let A be any Hermitian operator and Ajj; (A) be 
the eigenvalues of A arranged in decreasing order. From 
Cauchy's interlacing theorem [l[ (CIT) we have 

At(npn) < \f(p). (6) 

The condition tr[pLT] = ^i(p) implies 

d d 

^At(n P n) = ^Ai( P ) (7) 

i=l i=l 

which can only be satisfied if \\(JlpT[) = \\{p) for all 
i = 1 , . . . , d. This proves 

/ Ai(p) \ 

Pn =UpU= 4(P) (8) 



and p has the block structure 

Since ^jifLl = Pjji we have that p must be 

block-diagonal [2j; hence: 

»=( m n P n)- <»> 

□ 



Proofs of theorems 

[Theorem] Necessary and sufficient conditions for a 
Heisenberg Hamiltonian H to be mapped into a prototype 
model by a local gauge transformation are 

2a) All coupling vectors Jij, have the smallest absolute 
value component along the same axis; 

2b) For any pair of spins a and b, ir(p(a — > b)) is inde- 
pendent of the path p from a to b. 

When the second condition is met, the dependence 
of 7r on the path can be dropped. Moreover, we 
have 7r(a — > c) = 7r(a — > b)Tr(b — > c). The conditions are 
clearly satisfied for the prototype model. The subgroup 
H of G which preserves the Heisenberg structure of the 
Hamiltonian when acting on individual sites is gener- 
ated by the operations P^(-) = P(S^ ■ S^ 1 ) which all 
have diagonal 0(3) representations. Notice that each 
P M combines of an SU(2) transformation and a parity 
inversion. The representation of P 11 is given by i? M = 
diag((— l) 5 ^, (— l) 5 "^, (— l) 5 '^). Under these transfor- 
mations the component in J with smallest absolute value 
cannot be changed without H losing the Heisenberg form. 
Then, the path-independence of 7r(a — > b) is clear from 
the fact that if P M acts either on a or on b it will flip 
the sign of all paths from a to b. On the other hand, 
if P M acts on any other site, it will change none. This 
shows that 2a-b) are necessary conditions. To prove that 
they are sufficient, start from any spin a and apply to 
any other spin b the operation Rb = ir(a — > b). Then, 
for any pair of interacting spins i,j, the resulting new 
couplings are given by J(j = RjJijRj = \Jij\, where 
\Jij\ = diag(|J? |, \Jfj\, \Jij\)- This ensures that the axis 
of minimum coupling is preserved under such transforma- 
tions, and that all interactions are ferromagnetic. This 
proves the theorem. □ 

Testing conditions (2a-b) is particularly simplified if 
we consider translationally invariant systems, which play 
a fundamental role in statistical physics. In this case the 
system is fully described by the elementary cell, and we 
can state the following. 



6 




FIG. 1: Two paths pi and P2 going from a to b, passing 
through ci [solid] and C2 [dashed] respectively. 

A translationally invariant system obeys conditions (2a- 
b) if and only if they are satisfied for the elementary cell. 

The only if implication follows naturally. For the if im- 
plication consider two adjacent elementary cells, C\ and 
C2 [Fig- HI- Condition (2a) is automatically satisfied. To 
see that (2b) is also satisfied let a £ C\ and b £ C 2 , and 
consider two paths p\ , p 2 where pi (a — > b) goes from a to 
b passing through an intermediate site Ci £ C\ n C 2 . Then 

7t(pi(a 6)) =7r(a ci)[7r(ci c 2 )] 2 7r(ci 6) 
= 7r(p 2 (a -> 6)). 

The proof is concluded by realizing that this argument 
trivially extends to other elementary cells. □ 
Remarkably, a generalized form of Toulouse criterion 
for testing conditions (2a-b) on the elementary cell can 
be derived. 

If a model satisfies 

tt(£) = +1. (11) 

for every loop £ in the elementary cell and all coupling 
vectors have the smallest component along the same axis, 
then it is an INES model. 

If conditions 2a-b) are satisfied then the model is an 
INES model. For translationally invariant systems, it suf- 
fices that the conditions hold within the elementary cell. 
Hence, it only remains to prove that Eq. (6) is necessary 
and sufficient for the elementary cell to obey 2a-b). 

Any two paths px(a — > 6), P2(a — > b) can be merged 
into a closed loop I = p\{a — > b) U p 2 (a —> b), hence 2b) 
implies ir(l) = [ir(a —> b)] 2 = +1. On the other hand, 
if any closed loop yields ir(£) = +1, then any two paths 
Pi(a — > b), p 2 (a — > b) must have 7r(pi) = 7r(p 2 ). D 

Numerical evidence of the sufficient conditions 
(Conjectures 1 and 2) 

In order to test Conjectures 1 & 2 we have generated 
several Hamiltonians without geometric frustration, and 
checked that they satisfy the equality / = E in Eq.(3) 
in the main text for all pairs of spins. Instances of the 
prototype model have been generated, where the gauge 



was fixed for the unitary sector of the gauge group, since 
this symmetry is trivially proven analytically. In order 
to check Conjecture 2, a random gauge transformation 
consisting of pure parity transformations was applied. 
In total 278400 geometric frustration-free Hamiltonians 
have been checked by exact diagonalization. In all those 
without global ground-state degeneracy, no instance was 
found which failed to satisfy f = E for all pairs of spins. 

In the following we provide more details of our pro- 
cedure. We have evaluated 278400 spin-1/2 models on 
random graphs with n = 4, 6 and 8 spins (sites). All 
the numerical tests have been done using the program 
MATHEMATICA 6. We have generated random models 
of the XYZ type (general anisotropic Hciscnbcrg), XXZ 
and XXX (Hciscnbcrg) satisfying conditions 2a-b), both 
with homogeneous and inhomogeneous couplings. The 
models have been generated according to the follow- 
ing procedure. First, a connectivity graph (V,£) with 
|V| = n vertices connected by edges £ is generated. For 
each edge (ij) £ £ , an interaction of the types XYZ, 
XXZ or XXX was generated. Depending on whether 
an homogeneous or inhomogeneous model was being gen- 
erated, the interactions were obtained independently, or 
the same interaction was generated for all edges respec- 
tively. Finally, the identity or a parity inversion was ap- 
plied to each spin, randomly with probability 1/2. At 
each value of the edge probability q, and for each type 
of model, the method exhaustively evaluates all possi- 
ble connected simple graphs with n vertices. Specifi- 
cally, given a simple connected graph (V, £ ), the resulting 
Hamiltonians are, for the homogeneous case 

H hom = - ^ E f-) Ii8liS i s i' ( 12 ) 

v£{x,y,z] (ij)e£ 

and for the inhomogeneous case, 

#inhom = - {-) Xi&Xi JjjSj S" , (13) 

ve{x,y,z} (%j)es 

where the equiprobable binary random variables 

and each interaction was obtained by generating a ferro- 
magnetic coupling vector with Jf- , JV- and J?- , each 
component according to a flat p.d.f. (probability distri- 
bution function) \dp{J) = dJ] between and 1 for the 
general XYZ case. Ensuring condition 2a) was achieved 
by placing the smallest component on the z axis, and 
leaving the other two components in their original order. 
Notice that the choice of the z axis for the smallest com- 
ponent is arbitrary and the two other choices are equiva- 
lent up to a global SU(2) transformation. For the XXZ 
model the J x and J v components were symmetrized by 
replacing the obtained values by their averages. Finally, 
the XXX interactions were obtained by simply generat- 
ing random numbers between and 1. 

The connectivity graphs were obtained as follows. For 
n — 4 and n = 6, we tested all connected simple graphs 
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with n vertices. For n = 8, and due to the very large 
number of inequivalent connected simple graphs (11117), 
we iterated the Erdos- Renyi model Q to generate simple 
connected random graphs. In order to test all different 
kinds of connectivity (ranging from systems with few in- 
teractions such as the one-dimensional spin chain up to 
the fully connected model), the edge probability q was 
run from 0.1 to 1 with increments of 0.1 steps. At each 
value of q, and for each type of model (XYZ, XX Z or 
XXX), 2000 instances were generated, totaling at 18000 
models for each type of interaction. 

In Tables HI HT1 and [TTT1 we present our results, showing 
the total number of models generated for each type of 
interaction considered. The number of accepted Hamil- 
tonians (those for which the condition / = E is satisfied 
for all bonds (ij) £ £), rejected Hamiltonians (those vi- 
olating / = E for some bond) and degenerate Hamilto- 
nians (those with globally degenerate ground state) are 



shown for each interaction class. As shown by the ta- 
bles, no Hamiltonians have been Rejected, meaning that 
no instance has been found violating conjectures 1 & 2. 
This finding lends strong support to the general validity 
of the conjectures. 
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n # graphs Accepted Rejected Degenerate 



4 


6 


6000 








6 


112 


22400 









8 11117 18000 



4 6 6000 

6 112 22400 

8 11117 18000 



TABLE I: Results for XYZ models with homogeneous 
(above) and inhomogeneous (below) interactions. 



n # graphs Accepted Rejected Degenerate 



4 


6 


6000 








6 


112 


22400 








8 


11117 


18000 









4 6 6000 

6 112 22400 

8 11117 18000 



TABLE II: Results for XXZ models with homogeneous 
(above) and inhomogeneous (below) interactions. 



n # graphs Accepted Rejected Degenerate 



4 


6 


2247 





3753 


6 


112 


7072 





15328 


8 


11117 


4902 





13098 



4 


6 


2262 





3738 


6 


112 


7088 





15312 


8 


11117 


4908 





13092 



TABLE III: Results for XXX models with homogeneous 
(above) and inhomogeneous (below) interactions. 



